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Technical Report on the Physics of Heavy Ions (1989-1990) 

During the last year we made significant progress on understanding 
the enhancement of 0 + outflow in the polar wind. We completed studies 
on polar wind ion heating due to kinetic ion beam instabilities and the 
effects of such ion heating on the outflow of 0 + in the polar wind. 

DE 1 observations of 0 + and H + ion beams (~ 10 — 100 eV) in the 
polar cap at altitudes of ~ 2.5 — 3.5 R e found in a previous joint study 
with NASA Marshall Space Flight Center and Lockheed Palo Alto [Chen 
ei al., 1990] motivated the kinetic study of ion beam instabilities. The 
upflowing field-aligned 0 + and H" 1 " beams which may have been formed 
as upflowing-energetic ions from the dayside cusp E x B convect into 
the polar cap provide free energy to the polar cap plasma (see figure 1). 

We first examined the linear instabilities associated with an 0 + and H+ 
polar wind plasma in the presence of 0 + and H + beams for a range of 
0 + /H + beam densities, T e /T{, and ion beam speeds. Then, we studied 
the nonlinear heating of the polar wind ions using numerical simulations. ' 

We modelled the 0 + and H + polar wind ions by isotropic Maxwellian , 
distributions and the electrons, 0 + beams and H"*" beam ions by drifting 
Maxwellian distributions. The electron distribution drifted so that there 
was no zeroth-order current in the plasma. 

It was found that in the presence of 0 + and H + beams, the slow 
ion acoustic and slow ion cyclotron modes can couple to the normal 
modes of the plasma (for example, the 0 + and H + cyclotron modes 
for perpendicular propagation) and be driven unstable. When the 0 + 
beam density is dominant over the H + beam density, which is sometimes 
observed [Chen et al ., 1990], the slow 0 + acoustic and/or slow 0 + 
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cyclotron modes are triggered unstable while the H + beam modes are 
damped. Figure 2 shows the dispersion relation, frequency normalized 
by the 0 + gyrofrequency versus wavenumber normalized by the plasma 
frequency over the speed of light, of both the real part of the frequency 
and the growth rate for such a case. Note that the slow 0 + beam modes 
are unstable for a large range of wavenumbers due to their coupling 
with several O**" and H+ cyclotron harmonics associated with the polar 
wind. As T e approaches T{, the maximum growth rate of the slow ion 
acoustic mode decreases due to damping of the waves by the ions. At 
low beam speeds, the growth rate due to the coupling of the slow 0 + 
beam modes with the 0 + cyclotron harmonics are dominant while at 
larger ion beam speeds, the growth rates of the coupling with the H + 
cyclotron harmonics are dominant. 

Particle simulations showed that when the 0 + beam density is dom- 
inant over the H + beam density, the slow 0 + acoustic mode predicted 
by linear theory to be unstable is excited and can resonantly heat both 
the 0 + and H + polar wind ions in the perpendicular direction. In fact, 
the 0 + polar wind ions are preferentially heated over the H + polar wind 
ions when the O" 1 " beam density is dominant over the H^" beam density. 
Figure 3 shows the time histories of the normalized perpendicular H + 
and 0 + temperatures and electrostatic energy. The polar wind ions are 
significantly heated at the expense of the electrostatic wave energy. 

This ion heating could have very important implications on the po- 
lar wind outflow if it can alter the polar wind ions’ scale heights. We 
examined the effects of the kinetic ion heating on the outflow of the po- 
lar wind ions from the ionosphere using a time-dependent hydrodynamic 
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model- We developed a numerical code to solve the 0 + and H + continu- 
ity and momentum equations in a flux tube from ionospheric (200 km) 
to magnetospheric altitudes (~ 6 - 7 Re)- We included the effect of ion 
heating by allowing for the altitudinal variation of the ion temperatures 
in the momentum equation. The ion temperature profiles were specified 
based on the ion heating characteristics found from our previous kinetic 
simulations. We assumed that the heating occurred above 1500 km and 
increased to a saturated value of temperature that we obtained directly 
from the kinetic simulation study (see Figure 3). Since electrons were 
not significantly heated due to the ion beam instabilities, the electron 
temperature was kept constant. 

First, we studied the characteristics of the dynamical polar wind 
without ion heating. We simulated a flux tube on closed field lines 
(u t - = 0) that suddenly became open at t = 0. Figure 4 shows the density 
and velocity profiles of the H" 1 " and 0+ polar wind ions at various times 
in the absence of ion heating. In comparison to the H + ions, there is not 
much outflow of 0 + in the polar magnetosphere. Although there may 
be temporal enhancements of the 0 + fluxes, they last for only about 
15 min. and are due mainly to the initially large upward gravitational 
force when the flux tube suddenly becomes open. Clearly, there is not 
much sustained escape of 0 + ions from the polar ionosphere without 
ion heating. 

We then included the effects of ion heating. To gain some physi- 
cal insight, two limiting cases were considered: preferential heating 
(Tjj > To) and preferential 0 + heating ( Tq > T H ). These two cases 
corresponded to when the 0 + beam density is comparable to and domi- 
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nant over the H + beam density in the kinetic study, respectively. Figure 
5 shows the steady state results for these cases superimposed on the no 
heating case ( Tq = Tg). When the H + ions are heated, the H + velocity 
increases. This is because the heating increases the upward pressure 
gradient force and accelerates the H"*” ions at higher altitudes. However, 
the density does not alter much because the amount of H + which 
can escape the polar ionosphere is limited by charge exchange of neutral 
H with 0+ in the collisional ionosphere. Since ion heating occurs where 
the plasma is collisionless, H + heating does not affect the outflow that 
much For the 0^ ions, there is not only an increase of O'*" velocity but 
a dramatic increase of 0 + scale height above the heating region. This is 
because the 0 + ions were essentially gravitationally bound and the 0 + 
heating which increases the upward pressure gradient can help enhance 
the 0 + outflow. It does not matter whether the 0 + or H + ions are 
preferentially heated. As long as there is significant 0 + heating, the 
0 + escape from the polar ionosphere is enhanced. Thus, we showed 
how 0 + heating can lead to enhanced polar wind 0+ fluxes in the polar 
magnetosphere. 
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Fig. 1. A shematic of the polar magnetosphere which shows 
the polar wind flowing out of the polar ionosphere and the 
direction of the E x B convection of upflowing energetic ions 
from the dayside polar cusp region into the polar cap region. 
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Fig. 2. The dispersion relation when the 0 + beam is dominant 
over the H + beam. The modes and their growth rates are 
distinguished by different legends. (SOA = slow 0 + acoustic, 
SOC(n) = nth slow O” 1 " cyclotron harmonic). 
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Fig. 3. The time history of the perpendicular (a) H ' and 
(6) 0 + temperatures normalized by their initial perpendicular 
temperatures. The time history of the electrostatic energy is 
shown in (c). 
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Fig. 4. The altitude (km) versus (a) H"*~ density, ( b ) H~*" Mach 
number, (c) 0 + density, ( d ) 0 + Mach number for later times 

in the simulation when both H + and 0 + ions are included in 
the plasma. 
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Fig. 5. The (a) H + density, (6) H + Mach number, (c) 0 + den- 
sity, (d) 0 + Mach number profiles for 3 cases: no ion heating, 
preferential 0 + heating, and preferential H + heating. 
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